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ABSTRACT 

We investigate the brightness distribution expected for thermonuclear explosions that might 
result from the ignition of a detonation during the violent merger of white dwarf (WD) bina- 
ries. Violent WD mergers are a subclass of the canonical double degenerate scenario where 
two carbon-oxygen (CO) WDs merge when the larger WD fills its Roche-lobe. Determining 
their brightness distribution is critical for evaluating whether such an explosion model could 
be responsible for a significant fraction of the observed population of Type la supernovae 
(SNe la). We argue that the brightness of an explosion realised via the violent merger model 
is mainly determined by the mass of 56 Ni produced in the detonation of the primary CO WD. 
To quantify this link, we use a set of sub-Chandrasekhar mass WD detonation models to derive 
a relationship between primary WD mass (towd) and expected peak bolometric brightness 
(Mboi). We use this mwD--Wboi relationship to convert the masses of merging primary WDs 
from binary population models to a predicted distribution of explosion brightness. We also in- 
vestigate the sensitivity of our results to assumptions about the conditions required to realise a 
detonation during violent mergers of WDs. We find a striking similarity between the shape of 
our theoretical peak-magnitude distribution and that observed for SNe la: our model produces 
a Mboi distribution that roughly covers the range and matches the shape of the one observed 
for SNe la. However, this agreement hinges on a particular phase of mass accretion during 
binary evolution: the primary WD gains ~0.15 — 0.35 M from a slightly-evolved helium 
star companion. In our standard binary evolution model, such an accretion phase is predicted 
to occur for about 43% of all binary systems that ultimately give rise to binary CO WD merg- 
ers. We also find that with high probability, violent WD mergers involving the most massive 
primaries (>1.3M Q , which should produce bright SNe) have delay times <^500 Myr. 

Key words: hydrodynamics - radiative transfer - methods: numerical - binaries: close - 
supernovae: general - white dwarfs 
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1 INTRODUCTION 

It is widely agreed that the progenitors of Type la supernovae 
(SNe la) are thermonuclear explosions of carbon-oxygen (CO) 
white dwarfs (WDs) in binary systems. However, it is still not 
clear whether their companion stars (secondaries) are hydrogen- 
rich (main sequence or giant stars, typically called 'single degener- 
ate' scenarios), helium-rich (helium stars or helium-rich WDs) or 



* E-mail: a j r@mpa-garching . mpg . de 



other CO WDs which merge with the primaries (typically called 
'double degenerate' scenarios). 

I Sternberg et al.l d201ll) analysed high-resolution spectra of 
35 SNe la, finding that more than half of the spectra exhibit a 
blueshifted Na I D absorption feature. The authors attribute this 
feature to circumstellar material (blown off of the secondary prior 
to the SN explosion), suggesting that at least 20% of SNe la 
in s piral galaxies orig inate from single degenerate scenarios (see 
also IPatat et all 1201 ll) . On the other hand, several recent stud- 
ies have highlighted challenges to single degenerate scenarios 
and have instead tended to favour double degenerate models. 
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For example, the lack of emission associated w i th the super- 
nova overrunning the companion star (Kasen 2o Tol: iHayden et all 
l20ld:lGaneshalingam et alj|201 lLlBianco et alj|20ilh. the absence 
of radio detection s dHancock et alj 1201 ll ; Chomiuk et al .1 1201 2t 
lHoreshetalJl2012h and the lack of unambiguous identifications 
of companion stars in supern ova remnants dRuiz-Lapuente et all 
|2004 iKerzendorf et al.|[20Qgh argue against models in which the 
comp anion is a red giant star at explosion (see also iBloom et al .1 
HoH where a single degenerate scenario origin for SN 201 lfe is 
likely ruled out). In addition, a number of recent studies find that the 
delay time distribution (DTD) of SNe la follows a power-law shape 
t s , where -16 < s £ -1 jMaoz et alj|201(t iGraur et aljkfjl it 
iBarbarv et ai1l2012l : ISand et alj l2012h ." Such a distribution is ex- 
pected if the progenitor population of SNe la is dominated by dou- 
ble degenerates. These studies, among others, have fuelled specu- 
lation that there must be at least two different progenitor scenar- 
ios leading to SNe la. However, the question of which channel(s) 
dominate(s) the pro duction of 'normal' SNe la l lBenetti et al.l200 3l 
iBranch et alj2009l) and determine(s) the observed SN la brightness 

distribution remains open. 

In a series of recent papers jPakmor et al1l2010l . 1201 ll. |2012| ; 

iRopkeetal .1120121) . a new explosion model for SNe la produced 
by double degenerate (CO+CO WDs) mergers was investigated. In 
this violent merger model, a prompt detonation is triggered dur- 
ing the merger itself, leading to a thermonuclear explosion of the 
system. Although it is not yet established which parameters of 
the merging system are most crucia l for the formation of such a 
detonation (see e.g.. |Pan et alj|2012h , synthetic spectra and light- 
curves computed for several realisations of the violent merger sce- 
nario are a good match to maximum light observations of SNe la 
jPakmor et al]|2010ll2012h . The level of agreement with observa- 
tions is comparable to that found for other scenarios jRopke et all 
l2012h . 

In contrast to Chandrasekhar mass explosion models, where 
the peak brightness is driven by stochastic processes connected to 
the formation of deflagration and detonation flames o r other prop- 
erties of the exploding Chandrasekhar mass WD ( e.g. lKasen et alj 
120091 : iBlondin et al.l [20TTI ; ISeitenzahl et alj|201ll) . in the violent 
merger scenario there is a direct correlation between fundamental 
parameters of the progenitor system and the SN luminosity. Specif- 
ically, to first order, the violent merger model can be interpreted as 
a mechanism to ignite a detonation of a su b-Chandrasekhar- mass 
WD (the primary CO WD) as described in lSim et al.l d20ld) (see 
also lFink et al.ll201oh . Although the merging secondary WD is also 
consumed by the thermonuclear flame of the detonation, it does not 
produce any iron group elements when it is burned. Thus, the ab- 
solute brightness of a SN la explosion depends on the mass of the 
primary CO WD. 

The brightness distribution of SNe la can be observed directly, 
with a typical peak brigh tness of ~ — 19 mag ^Richardson et alj 
l2002l :l Benetti et alj l2005t) and a spread in brightness of ~1 mag. 
Recently, Li et alj d201 ll) presented the brightness distribution of a 



volume-limited sample of SNe la. Any theoretical scenario explain- 
ing SNe la that claims to account for a large fraction of observed 
events must be able to reproduce the observed brightness distribu- 
tion. The goal of this paper is to quantify the brightness distribution 
of SNe la from the violent merger scenario. 

Our study requires three main inputs. First, we need predic- 
tions of the distribution of properties of merging CO+CO WD 
pairs. These we obtain from previously published StarTrack 
(Belczvnski et al ]|2002l . l2008l) binary population synthesis models, 
which are described in Section|2] Second, we must consider which 



characteristics of a merging system will be required for the explo- 
sion mechanism to work. This specification will be made based 
on the mass ratio of the merging system, as described in Sec- 
tion [3] Third, we need to estimate the brightness of an explosion 
based on the properties of the merging WDs. This is discussed 
in Section [4] In Section [5] we present our results based on the 
standard model us ed for the double degenerate scenario (DDS) in 
iRuiter et"al] d201 lh and comment on the comparison with observa- 
tions. We discuss important implications and uncertainties of our 
modelling in Section|6]before drawing conclusions in Section|7] 



2 BINARY POPULATION SYNTHESIS 

2.1 Mass exchange for StarTrack binaries 

We consider the p opulation synthesis calculation results from the 
standard model of IRuiter et al. I J20T l|), whose data were obtained 
with the StarTrack population synthesis binary evolution code. 
In this subsection, we briefly explain how stellar mass loss/gain is 
handled over the course of binary evolution in StarTrack. 

Unstable mass transfer. Outside densely populated environ- 
ments - such as the cores of globular clusters - binary stars 
that form close doubl e WDs must underg o at least one common 
envelope (CE) phase dlben & Livid 1 19931) , in which the size of 
the binary orbit decreases (often significantly) upon expulsion of 
the mass-losing star's envelope. Currently, CE ejection remains a 
poorly understood process in astrophysics. Observationally, it is 
difficult to study since it is a short-lived evolutionary phasj]] and 
detailed theoretical modelling is extremely challenging given the 
large spatial and temporal scales that must be taken into account to 
properly approach the problem jRicker & Taarr]|2012l) . However, 
despite our collective ignorance, the CE phase cannot be ignored 
in binary evolution modelling, and so a simplified parametrization 
of the CE phase is often employed in binary population synthesis 
codes. 

In the standard model of IRuiter et af] d201 lh . the CE 
parametrization o?ceA = 1 was use d. In this CE pa rametrization, 
often coined the 'energy formalism ' dWebbinklll984l) . it is assumed 
that the binary can expel the envelope without merging if the bind- 
ing energy of the mass-losing star's envelope is not greater than 
the orbital energy of the binary wit h some efficiency factor oce ; 
A-Ebin ^ Qc e Ai7 or b (but see also llvanova & Chaichenetsll201ll ; 
llvanovallioTl for an overview and recent work on t he CE prob- 
lem). S ince our calculations are based on the results of lRuiter et all 
J20I l|), we assume as they do that the transfer of energy from the 
orbit to the removal of the envelope during CE is fully efficient, 
and thus q?ce = 1. The binding energy is inversely proportional 
to the term A, which is a parametr ization of the structure of the 
mass- losing star (e.g. lde Koollll990|) . A is somewhat uncertain, but 
is often assumed to be on the order of unity in population synthe- 
sis calculations for low- and inter mediate-mass stars on the main 
sequence (MS) or the giant branch dDewi & Taurill2000l). We take 
the CE parametrization with qqe = 1 and A = 1 from Ruiter et all 
j201 ll) for the current study. 

Stable mass transfer. For Roche-lobe overflow (RLOF) on 
to non-degenerate accretors, mass transfer is assumed to be non- 
conservative. An efficiency parameter - which characterizes the 
fraction of matter transferred to the companion - is adopted (termed 



1 Though the sample of post-CE binaries (PCEBs) is becoming larger 
lZorotovicetalj201lh . 
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Figure 1. The mass distribution of the primary WD for all CO+CO WD 
binaries at the time of the merger from the Ruiter et al. (2011) standard 
calculation (solid histogram). The dashed histogram shows the mass distri- 
bution of the primary when it becomes a CO WD (i.e. neglecting any mass 
accretion which may occur after the primary WD has formed, see text). 
The numbers on the y— axis are not scaled, and show the actual numbers 
of CO+CO WD mergers produced in a Hubble time in the simulation of 
4.5 X 10 6 zero-age main sequence binary stars. 



in Equation (33) of iBelczvnski et al] l2008h . In iRuiter et al] 



20iTh. f. 



0.5 was adop ted for the standard calculations (see 
also Mennekens et aljfeoiol , where the same parameter is termed 
'/?'). Thus only 50% of the matter lost by the donor is gained by 
the accretor (Eddington-limited) while the rest is lost from the bi- 
nary system with the specific angular momentum of the orbit (see 
Appendix A for more discussion). 

For WD accretors, the StarTrack code takes into account a 
number of published resu lts on mass accretion on to CO WDs for 
hydro gen-rich accretion jPrialnik&Kovetj|l995 " 



9ll; iNomoto et al] 

l2007hrand helium-rich accretion dKato & Hachisulll999ll2004 see 



also lTaamlll980l ; lGarcia-Senz et al.ll 19990 . In this case, the fraction 
of mass gained by the accretor is dependent upon the rate at which 
mass is transferred from the donor and the mass of the accreting 
WD. Mass that is lost from the binary carries with it the specific 
angular mo mentum of the accretor. A more detailed description can 
be found in lBelczvnski et al] ( 120081 section 5.7.2). 



2.3 Binary evolution and example 

To understand the origin of the double-peaked mwD -distribution, 
we have analysed various binary evolution histories for representa- 
tive systems in the StarTrack simulations. The formation of the 
secondary bump at ~ 1.1 Mq in the solid histogram of Fig. 1 (and 
more generally the larger average primary WD masses) is associ- 
ated with CO WD primaries that accreted ma terial from a slight ly 
evolved naked helium star companion (k=8 in iHurlev et alfeOOOh . 

To be more specific, we present a typical evolutionary chan- 
nel along with the relevant physical processes that may potentially 
lead to the formation of a CO+CO WD merger that contributes to 
the bump at ~1.1 Mq in Fig.Q] In our example, the initially more 
massive star on the zero-age main sequence (ZAMS) will become 
the first formed WD and the more massive WD at the time of the 
merger, though this is not always the case for all evolutionary chan- 
nels. However, throughout the paper, we define the primary (with 
mass m p ) as the star which ends up as the more massive WD at 
the time of the merger, and the companion is the secondary (with 
mass m s ). We illustrate the onset of the most important evolution- 
ary stages in Fig. [2] with Roman numerals (I-X), where time t, 
masses m p and m s , and separation a are indicated in Myr, Mq and 
R©, respectively. We also show the mass evolution of the binary in 
Fig. [3] 

Our example begins with two MS stars on a circular orbit with 
an orbital period of 8.3 h (a=37 Rq, stage I). On the ZAMS, m p 
and m s are 5.65 and 4.32 Mq, and their corresponding radii are 
2.84 and 2.43 Rq, respectively. At 79 Myr the primary evolves 
off the MS and fills its Roche-lobe, leading to RLOF while in the 
Hertzsprung Gap (stage 11)0 The mass ratio is fairly close to unity, 
so RLOF is stable (/ a = 0.5) though it initially proceeds on a 
thermal time-scale. Over the course of RLOF, which lasts 0. 1 5 Myr, 
there is a mass ratio reversal, where the primary becomes the less 
massive star. The mass transfer (from less massive to more massive 
star) and expansion of the primary drive an increase in separation to 
233 Rq , the primary shrinks back inside its Roche-lobe and RLOF 
ceases. By the time RLOF has stopped, m p = 1.02 Mq while 
m s has increased to 6.62 Mq. The primary continues to evolve up 
the red giant branch, though by this point most of its hydrogen- 
rich layers have been removed in RLOF. Once the hydrogen shell- 
burning phase is complete, the primary becomes a naked helium 
star (core of red giant). At t=99 Myr, the primary becomes a slightly 
evolved naked helium star and increases significantly in size. At 
t=102 Myr, the primary fills its Roche-lobe and a second phase of 
stable RLOF ensues (stage III). Helium-rich material is transferred 
to the secondary on a nuclear time-scale for 0. 14 Myr (as before 
/ a = 0.5). Mass transfer ceases at a separation of a = 318 Rq, 



2.2 Primary mass distribution 

Fig. H shows the distribution of primary (more massive WD) 
ma sses for merging CO +CO WD pairs at the time of merger from 
the IRuiter et al] ( 1201 ll) «ceA = 1 calculation (solid histogram). 
Note that this mwD distribution has two peaks: a prominent peak 
at ~0.65 Mq and a smaller bump at ~1.1 Mq. For comparison, 
we show the distribution of primary masses for the same popula- 
tion of merging CO+CO WD pairs recorded at the time of primary 
WD birth (dashed histogram); this distribution has only a single 
peak at ~0.65 Mq. This 'mass shift' to larger masses is the result 
of a phase of mass accretion when the donor star is in the helium- 
burning phase, which we discuss in detail in Section [2731 As we 
will discuss in Section [4] this has important implications for our 
predicted brightness distribution of SNe la. 



evolution, the primary has converted much of its helium core into 
CO, and upon completion of helium burning at t — 102 Myr, it 
becomes a CO WD (stage IV). 

At t=l 15 Myr, the primary begins its journey up the red giant 
branch. Once the secondary reaches a radius of 123 Rq (orbital pe- 
riod of 140 d), a CE takes place (stages V, VI), leaving behind the 
CO WD and the (naked) core of the red giant, a decreased separa- 
tion (orbital period of 261 min) and m B — 1.27 Mq (stage VII). 
At t — 127 Myr, the secondary becomes a slightly evolved naked 
helium star. Soon afterwards it fills its Roche-lobe which leads to a 
third RLOF phase (stage VIII). 



2 See [Belczynski et al. 12008, sect. 5.6) for a description of how stellar 
rejuvenation is treated in STARTRACK. 
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The mass transfer rate initially proceeds on a thermal time- 
scale and is fairly high (Mtran = 8.8 x 1CT 5 M yr -1 ). For the 
given WD accretor mass, such a mass transfer rate is above the 
critical threshold that leads to 100% burning efficiency (the criti- 
cal v alue is 4.57 x 10~ 7 M q yr -1 for 0.84 M WD, interpolated 
from lKato & Hac hisu 2004, Equation 6). However, the mass trans- 
fer rate is super-Eddington, and so the accretion rate is capped at 
the Eddington accretion rate (M aC cr = 2.0 x 10 _5 M Q yr" 1 ; Ed- 
dington accretion lasting 0.002 Myr)0 The mass transfer rate soon 
after enters the nuclear time-scale regime, and as it decreases he- 
lium burning continues to be stable on the CO WD. Once the mass 
transfer rate drops below Mtran = 1-74 x 10~ 6 M yr -1 (still 
t = 128 Myr; WD mass 1.08 Mq), helium burning is no longer 
stable (Equation 4, Kato & Hachisu 2004) and thus not fully ef- 
ficient. The primary continues to accrete another 0.11 Mq (total 
mass transfer including Eddington-limited phase lasts 0.15 Myr), 
though RLOF shuts off when the orbit has expanded to a = 1.91 
Rq (orbital period of 5.2 hr). At this point, m p = 1.19 Mq and 
m s — 0.78 Mq, while the secondary's helium core is almost com- 
pletely burned to CO (m COIC = 0.76 Mq). Within 10 5 yr (after los- 
ing ~0.001 Mq in stellar winds) the secondary becomes a CO WD 
(stage IX). The two WDs evolve toward smaller orbital separations 
over the next Gyr under the influence of gravitational wave radia- 
tion, and merge at t = 1259 Myr (stage X), with a combined mass 
of 1.96 Mq. Such a delay time (1259 Myr) is typical of double de- 
generate mergers - they are predicted to produce more SNe la at 
early times, peaking at a few hundred Myr though sti ll producing 
a sign ificant number of events >1 — 3 Gyr (e.g. fig. l. lRuiter et al .1 
120091) . 

The system described above is typical of an evolutionary chan- 
nel in which the critical accretion phase onto a WD (hereafter stage 
VIII, cf. Fig. 2) occurs . For the evolution ary model under consid- 
eration (oceA = 1 of iRuiter et al]|201ll) . we find that 43% of all 
CO+CO WD merger progenitors undergo such an accretion phase 
during which the primary WD (e.g. the first-formed and more mas- 
sive WD) accretes an additional 0.001 — 0.45 Mq after its forma- 
tion (in the majority of cases it gains ~0.2 ± 0.1 Mq). In addition 
to these 43%, 9% of merger progenitors encounter the opposite sit- 
uation: the secondary WD (e.g. the first-formed but less-massive 
WD) undergoes helium accretion^ The remaining 48% of progen- 
itors do not encounter a stable mass accretion phase while one of 
the stars is a WD. 

Different physical properties of a binary at birth (masses, mass 
ratio and separation/eccentricity) influence the ensuing evolution. 
We find that for the (43%) group of primary interest, typical ZAMS 
masses are somewhat large, peaking at ~4 — 5 Mq (though the 
entire range is 2 — 7 Mq), and initial mass ratios qo(= more mas- 
sive star / less massive star) range ~1 — 2 with a peak at ~1. For 
cases where the WD-accretion phase occurs on the secondary WD 
(9%), ZAMS masses are in the range ~2 — 3.5 M and thus re- 
sult in lower-mass WDs (in addition, qo for these systems peaks 
at ~1.4). For merger progenitors which do not undergo the WD- 
accretion phase (48%), ZAMS masses are mostly ~2 — 3 M , and 
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Figure 2. Flowchart showing a typical evolutionary channel where the pri- 
mary WD gains mass via RLOF from its helium star companion, prior to 
the double CO WD merger. The duration of the RLOF phases and other 
evolutionary details are given in the text. 



3 If such binaries lead to SNe la, we estimate ~few to 10 systems to be 
in this high-accretion phase and thus brig ht and potentially visibl e in the 
Galaxy (Galactic SN la rate of 0.007 vr~ 1 . iBadenes & MaodioLi section 
4). This number is within observational constraints, and in fact our example 
system exhibits physical properties similar to those of the helium nova V445 
Puppis iKato et alj200j) . 

4 These secondaries only gain a small amount of mass and they never be- 
come the more massive (exploding) WD. 



qo strongly peaks at unity. Hence, a large fraction of these stars 
evolve at roughly equal rates and in many cases, they reach the 
short-lived helium star phase at roughly equal times, so stage VIII 
does not occur. 
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Table 1. The peak bolometric magnitude derived from our ID sub- 
Chandrasekhar mass pure detonation models for CO WDs for a range of 
WD masses (mwD)- 
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Figure 3. Mass evolution for the binary system shown in Fig. [2] We do not 
show the last ~Gyr of evolution since the WD masses remain unchanged 
during that time. The common envelope event is marked by the drastic de- 
crease in secondary mass at 115 Myr, while the two vertical lines indicate 
the formation times of the CO WDs. 

3 CONDITION FOR PROMPT DETONATION 

Binary population synthesis calculations provide a predicted dis- 
tribution of physical parameters for merging CO W D pairs. How- 
ever, the prompt de tonation explosion mechanism of lPakmor et al.l 
J2010L 1201 lL l2012h is not expected to b e realised for all CO WD 
mergers. As discussed by IPakmor et al. I boid I20T l|), conditions 
favouring a prompt detonation are most plausible when the ratio 
of the secondary WD mass (m s ) to the primary WD mass (m p ) is 
high (approaching 1 ). In a study of mergers with m p = 0.9 Mq, 
IPakmor et alj d201 ll) determined that the critical mass ratio (q c = 
m B /m v ) above which explosion is plausible is q c ~0.8. The criti- 
cal q- value, however, may depend on the mass of the primary, likely 
having a value that decreases with increasing m p (for higher-mass 
primaries, the densities are higher and the merger is more violent 
making it easier to reach the con ditions needed for triggering a 
detonation, ISeitenzahl et alj|2009l) . Presently, sufficiently detailed 
simulations to determine how q c will depend on primary mass have 
not been carried out. Therefore, in this study, we will adopt a simple 
parametrization of q c : 

?c = min ( a8 (aSfe)"' L0) (1) 

where 77 is a parameter that we will use to explore how q c influences 
our findings. 



4 PRIMARY MASS/PEAK BRIGHTNESS RELATION 

Ideally, synthetic observables for a suite of three-dimensional WD 
merger models would be used to determine the relationship be- 
tween the system parameters and the explosion brightness. Unfor- 
tunately, such simulations are too computationally expensive for 
this to be practical in this study. However, as argued in Section Q] 
the peak brightness of a prompt detonation merger model is directly 
related to the mass of the more massive WD (m p ); in the merger, 
the 56 Ni is mostly synthesized in the detonation of the primary. We 



therefore estimate the brightness of a merger from a simple one- 
dimensional model for a detonation of a WD with mass equal to 
that of the primary. 

To quantify the relationship between the primary mass and the 
peak bolometric magnitude (Afboi), we carried out a series of sim- 
ulations for detonations of single CO WDs with a range of WD 
masses (mwD = m p ). These simulations were computed i n ex- 
actly the same manner as those described by ISim et aljj2O10l) . We 
constructed hydrostatic models for CO WDs and simulated deto- 
nations of these WDs using our SN la explosion code. In all cases, 
the detonations were centrally igni ted and their propagation fol- 
lowed using the level-set technique dGolombek & Niemeve r|[20o3 : 
iFink et al]|2010h . Nucleosynthesis in the explosions was computed 
in a post-processing st e p using our standard tra cer particle approach 
dTravaglio etai]|2004 ISeitenzahl et alj|2010h assuming an initial 
12 C/ 16 0/ 22 Ne composition of 47.5/50/2.5 % by mass (roughly 
appropriate for solar metallicity). For each explosion model we 
performed radiative transfer simulations with ARTIS Jsimll2007l : 
iKromer & Siml2009l) to pre dict synthetic light-cu rves (adopting the 
cd23_gf-5 atomic dataset of lKromer & Siml l2009h . 

Our chosen values of mwD and the derived values of Mboi are 
given in Table 1. Fig.|4]shows the mwD-Afboi relationship derived 
from the models. As expected, the models define a smo oth, mono- 
tonic relationship between brightness and mass (see e.g. lSim et al .1 
|2010|) . For the analysis below, we have made a polynomial inter- 
polation between the data points defined by the models to obtain a 
simple functional form for the mwD-Mb„i relationship. This inter- 
polated relationship is also shown in Fig. [4] (see Appendix B for 
coefficients of the fit). 

Recently, IPakmor et al.l d2012l) presented results from a com- 
plete simulation of a prompt-detonation model for the merger of 
a CO WD pair with m p = 1.1 and m s = 0.9 M©. The angle- 
averaged peak brightness for this model is indicated in Fig. [4] As 
expected, its peak brightness is very close (AAfboi ~ 0.2 mag) to 
that of our equivalent single WD detonation model. We also show 
th e angle-averag e d peak Mboi for the complete merger simulation 
of IPakmor etafl d2010h . A gain, the agreement with our relation 
is good, although the difference is slightly larger (AMboi ~ 0.4 
mag). These comparisons lend credence to the approach we adopt 
to estimate the peak brightness of WD merger models and suggest 
that our method is accurate to a level of several tenths of a magni- 
tude. 
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Figure 4. The mwD-^bol relationship derived from the pure detonation 
models in Table 1 . The rilled circles indicate the data derived from the mod- 
els and the solid line shows the polynomial interpolation we adopt in our 
analysis. The squares indicate the results from the full merger simulations 
described by Pakmor et al. (2010; subluminous) and Pakmor et al. (2012; 
M bol ~ - 19). 



Figure 5. The predicted distribution of peak bolometric magnitude (iWbol) 
obtained by comb ining the primary WD masses from the model of 
iRuiter et alj feOl ll) with our TnyvD-^bol relationship (Fig. [4}. The black 
line shows the results assuming that all mergers successfully lead to a 
prompt detonation. The coloured lines show the effect of imposing cuts 
on the (/-value of the merger (requiring that q > q c ) for different values 
of rj (see text). The grey histogram shows the absolute m agnitude distri- 
bution derived from a volume-limited sample of SNe la bv lLi et all 1201 lL 
see their fig. 5). To facilitate comparison of the shapes, the observational 
data have been scaled to match the number of systems in our calculation for 
T] = 0.75. 



5 RESULTS 

5.1 Distribution of brightness 

To convert the population synthesis data to a distribution of peak 
brightness, we first select all systems from the binary population 
synthesis calculation that lead to CO+CO mergers. From these we 
identify those systems that will explode based on our mass ratio cut 
(q > q c , and chosen values of 77; Section|3]l- Finally, we apply our 
mwD-Mboi relationship (Section |4) to derive the peak brightness 
for each system from its primary WD mass. We exclude from our 
analysis systems in which the primary WD is m p < 0.8 Mq, for 
two reasons. First, for such small primary masses the low densities 
make it extremely challenging to ignite the detonation. Second, at 
low mass the 56 Ni yield is expected to be too small to give rise to 
a bright transient even if a prompt detonation were to occur (see 
Table 1). This effectively removes the majority of systems asso- 
ciated with the low-mwD peak (~0.65 Mq) from our analysis: 
those systems are not promising candidates for the prompt deto- 
nation scenario (although it has been suggested that expl osion of 
suc h systems might still occur. Ivan Kerkwiik et al.|[2010i : but see 
also lShen et ai1l2012l) . 

Our computed brightness distribution is plotted in Fig. [5] for 
three representative values of the ^-parameter in our prescription 
for q c (see Equation 1). We have restricted 77 to be positive, ensur- 
ing that the critical q-value decreases with increasing m p (except 
for the case where r\ — 0, where q c = 0.8 for all values of m p ). We 
find that for 77^ 1.5, nearly all of the primaries at the high-mass end 
are included, so adopting larger values of r\ is unnecessary. Thus, 
we choose two values of 77 which bracket the physically-realistic 
extremes: 77 = and 77 = 1.5, and a third value 77 = 0.75 (the mid- 



point), simply for illustrative purposes^ We also show the result 
without any g-cuts (i.e. all population synthesis data). 

As shown in Fig. 5, our mwD-Mbd relationship converts the 
higher mass peak of the m p distribution (Fig. Q] solid line) into a 
brightness distribution with a clear peak around A/boi— — 19 mag, 
a short tail of brighter events, and a long tail stretching to fainter 
magnitudes. The amplitude of the peak and to some degree its lo- 
cation depend on our choice of 77. For 77 = 1.5 (or higher), the high 
Mboi part of the distribution is almost insensitive to the choice of 77 
but for smaller values (e.g. 77 = 0), a significant fraction of the high 
m p systems are excluded. However, 77 does not qualitatively affect 
our findings: the peak around Mt>oi~ — 19 mag is present in all 
cases - it is a consequence of the underlying shape of the m p dis- 
tribution in merging systems and the steepness of our mwD-Mboi 
relationship. 

5.2 Comparison to observations 

Recently, iLi et alj d201 ll) presented pseudo-observed luminosity 
functions for SNe la based on a sample of objects within 80 
Mpc from the Lick O bservatory Supernova Search (LOSS, see e.g. 
iLeaman et al. I l201ll . for details). We compare the shape of their 
volume-limited luminosity function to the results of our calcula- 
tions in Fig.[5] In this comparison, we include all 74 of their SNe la 
(regardless of sub-type and host galaxy type) and, for convenience, 
have arbitrarily rescaled the distribution to match the total number 
of SNe la to that of our calculation for 77 = 0.75. 

5 We note that for our evolutionary example of a violent merger formation 
channel described in Section 2.3 (^merger = m a /m p = 0.774/1.185 = 
0.653), the necessary criteria leading to a SN la are satisfied for 77 > 0.74: 
77 = 1.5 -> q c = 0.530; 77 = 0.75 -> q c = 0.651. 
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Figure 6. Primary WD mass (ni p ) vs. delay time up to 2500 Myr for Star- 
TRACK CO+CO WD mergers. The upper panel shows all systems while the 
lower panel only includes systems that pass the q c criterion with r\ = 1 . 5 
and m p > 0.8Mq. Note the different scaling for the two plots. It is evi- 
dent that the most massive primaries preferentially merge at prompt delay 
times: for both panels ^95% of systems with m p > 1.3 Mg have delay 
times < 1 Gyr. Mergers with m p > 1.3 M make up 0.2% of all CO+CO 
mergers and 0.5% of rj = 1.5, m p > 0.8Mq mergers. We also note that 
the first progenitors to explode are not from binaries with the most massive 
primaries (see text for details). 



This comparison suggests that the prompt detonation vio- 
lent merger scenario could give rise to a brightness distribution 
of SNe la that is fairly c ompatible w i th obs ervations. In particu- 
lar, the standard model of lRuiter et al.l d201 lh predicts a large frac- 
tion of explosions within a fairly narrow range of peak magnitude 
(AAfboi ~ 1 mag) centred around Afboi~ — 19 mag, as is observed. 

There are a number of important caveats to consider when in- 
terpreting the comparison in Fig. [5] First, the supernovae in the 
iLi etalJfcOl 1|) luminosity function have been corrected for Galac- 
tic extinction but not for host galaxy reddening, thus variations in 
extinction from object to object will be responsible for some com- 
ponent of the t ail to lower pea k magnitudes in the observed sam- 
ple. Second, the lLi et al I J20T H) data are not derived from complete 
sets of colour light-curves but are available only for the i?-band. 
However, this is not expected to significantly affect our results 
around the peak of our distribution : in th e complete merger sim- 
ulation presented by IPakmor et ail J2012h . the difference in peak 
bolometric and i?-band magnitudes was Afboi — Mr = 0. We 
note that A/boi — Mr is expected to become positive for sub- 
luminous mergers but for the systems of interest, this should not be 
large effect compar ed to the bin size in Fig.[5](in the simulation of 
IPakmor et alj2010L 

A/boi — Mr~0.4 mag.). Third, our analysis ne- 
glects any influence of observer orien tation on the explosion bright- 
ness - the IPakmor et all d2010L l2012h simulations suggest that ori- 
entation may cause variations of around ±0.4 mag. Nevertheless, 
although quantitatively important, none of the effects is expected to 
qualitatively change our result. 



Figure 7. STARTRACK smoothed DTD of all CO+CO mergers (green dots), 
CO+CO mergers where Af to t > 1.4M (blue solid), and CO+CO WD 
mergers which satisfy our criteria 77 = 1.5 and m p > 0.8Mq (magenta 
dashes). The small fluctuations are due to Monte Carlo noise. Over-plotted 
with red squares is the best-fit DTD derived from SDSS data, the three time 
bins are indicated by horiz ontal error bars: t < 0.42 Gyr, 0.42 < t < 2.4 
Gyr and t > 2.4 Gyr (see lMaoz et"aT]|2Q12L MMB12) The smooth black 
solid line shows the best-fit DTD from lGraur & Maoj 420121 GM12), with 
systematic uncertainties (thinner dotted lines). The right-axis shows units in 
SNuM: number of SN la per 10 10 Mq formed in stars per century. 

5.2.1 Brightness and delay time 

We note the important finding that, while fairly massive WDs with 
masses of ~1.2 M Q may explode at any delay time >150 Myr, 
the most massive ones (^1.3 Mq) only explode as prompt events, 
with rare exceptions (see Figure . This finding supports the ob- 
servational evidence for mo re luminous SNe l a occurring among 
younger stellar populations terandtetai1l2010h . At the same time 
however, the events with the shortest delay times (< 100 Myr) 
are not from the most massive primaries: these occur via rarer for- 
mation channels where (usually) the secondary initiates two CE 
events. These 'ultra-prompt' double CE systems do not undergo 
phase VIII mass accretion, though they are still massive enough to 
be considered SNe la progenitors in our model. The difference in 
evolutionary timescale between these double CE events and all of 
the other formation channels is the cause of the dearth of mergers 
at ~150 Myr in Figure(6] For the majority of the event s at later de- 
lay times only one CE is enc ountered, as is expected l lRuiter et al .1 
l2009l : lMennekens et alj2010h . 

5.2.2 Rates and DTD 

We have not attempted to calibrate the population synthesis event 
frequency of CO+CO violen t mergers to compare directly with the 
SN la rates o f lLi et aljf201lh . To do so would require us to adopt a 
model for the star-formation history associated with their volume- 
limited sample, which would introduce unwanted uncertainties. 
Nonetheless, we can compare our theoretical DTD to that from ob- 
servations. Merging CO+CO WD pairs are currently the favoured 
theoretical scenario for producing a relatively high number of SN 
la progenitors compared to other long-standing progenitor model s 
Jlben&TutukBvIl 19841 ; lYungelson etafll 19941 : iRuiter et alj|2009h . 
In Figure [JJ we show three theoretical DTDs for WD mergers as- 
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suming a binary fractiorlf] of 70%: All CO+CO WDs, mergers for 
which A/tot > 1.4 A/© (DDS), and systems which meet our crite- 
ria 77 = 1.5 and m p > 0.8 Mq (violent mergers). Additionally we 
show two of the most up-to-date obs ervationally-recovered DTDs 
from the literature: Maoz et al. I fcOlj SDSS II galaxies; fig. 1) and 
iGraur & Maozl d2012l . field rates; fig. 12). 

The DTD from violent mergers (magenta dashed line) com- 
pares rather well with the most recent observational DTD esti- 
mates. At long delay times (t > 2.4 Gyr) the vio lent merger 
DTD falls within the observational uncertainties of iMaoz et alj 
Also, the model a grees very well with the best-fit DTD 
from lbraur & Maozl d2012l) for delay times above 300 Myr, imply- 
ing that violent merg ers might be a domin ant contributor to SNe la 
in field galaxies (see lGraur & Maozll2012l . for a comparison of ob- 
servational DTDs from recent works). 

Up until now, there has been a notable discrepancy between 
observed rates of SNe la - in particular at early delay times - 
and the rates calculated from population synthesis. Typically, bi- 
nary population synthesis rates (as a function of delay time) were 
too low to explain the observed numbers. This discrepancy re- 
mains an op en issue, though some possible reasons h ave been dis- 
cussecQ (see lNelemans et al.l2012l ; rRuiter et alj|201ll . sect. 4. 1 ; see 
also Mennekens 2010 who point out that rotation may increase 
the likelihood for explosion in WD mergers). Howev er, the new 
observationally-derived DTD of IGraur & Maozl d2012l) has a lower 
amplitude (implying lower SN la rates) than those from previous 
observations. This brings our theoretical DTD, as well as some the- 
oretical DTDs from other grou ps, into better agreeme nt with obser- 
vations. It should be noted that lGraur & Maoa J2012l . sect. 6) men- 
tion that their best-fit power-law DTD might be underestimated, 
since it was derived using predominantly old galaxies. 



6 DISCUSSION 



We have found that, in the standard model of IRuiter et alj d201 lh . 
a certain formation channel is critical for creating massive primary 
WDs. In Section l2~3l we described the typical evolution of a binary 
which produces two merging CO WDs where the primary WD is 
particularly massive due to an important evolutionary phase: it ac- 
cretes (in that case 0.35 M©) from a helium-rich companion star. 
Within the model framework of lRuiter et alj ( l201ll) . such a phase of 
RLOF occurs for a significant fraction (43%) of the CO+CO WD 
merger progenitors. 

We have mapped the (exploding) primary WD masses for all 
CO+CO WD mergers on to magnitude- space using our mwD- 
Afboi relation, and have presented the results for a series of q- 
cuts (Fig. |5). Regardless of the g-cut used, the shape of our the- 
oretical A/boi distribution for merging primary WDs is quite strik- 
ing, sinc e it matches we ll with the observed Afboi distribution for 
SNe la o flLietal.lj20TTh . Our results indicate that double degener- 
ate mergers exploding in the violent merger scenario show excep- 
tional promise as progenitors of SNe la. 

However, despite these very encouraging results, we point out 
in the next two sub-sections some important issues to keep in mind, 



6 For stars of spectral type B3-B7, which are typical progenitors of vio- 
lent WD mergers, the bi nary fraction is estimated to lie between 40-100 % 
jHoffmeister et alj201 if). We do not consider higher multiplicity here. 

7 Observational samples which include massive elliptical galaxies tend to 
exhibit DTDs which have higher SN rates th an volume-limited samples, D. 
Maoz, private communication 2012; see also lMaoz et alj J2012I) . 



which could have a significant effect on the evolutionary outcomes 
for close binaries, in turn affecting the impact of these new findings. 

6.1 A key assumption: Common Envelope parametrization 

Different assumed CE parametrizations have a critical impact on 
which evolutionary channels lead to the formation of close bina- 
ries, as well as how often they form and how often they merge 
toominik et alj|2012l) . So, we now discuss how the results found 
above are sensitive to the choice of the adopted CE formalism and 
CE efficiency. 

Considering a lower CE efficiency such as qqe = 0.5 with 
A = 1 — > oceA = 0.5, which is often used in population syn- 
thesis, we find that the stage VIII accretion phase is indeed en- 
countered though not as often as in the standard model (33% vs. 
43%). Quantitatively speaking, the overall merger rates from pop- 
ulation synthesis calculations for this CE model are abou t a factor 
of two lower than the standard model ( IRuiter et al.ll2009l) . since a 
larger number of progenitors which make CO+CO WD mergers in 
the standard model will merge before a double WD is formed in 
a lower CE efficiency model. Qualitatively however, there is little 
difference between the two models, and in the qceA = 0.5 model 
we still expect such a 'stage VIII' accretion phase to occur leading 
(favourably) to more massive primary WDs. 

For t he case of very low CE efficiency, such as that inves- 
tigated in IRuiter et all J201 lh (a C E = 0.25 with A = 0.5 -> 
oceA = 0.125), stage VIII is almost never encountered: less 
than 1% of binaries go through this evolutionary channel (> 
99% of primary WDs do not experience accretion before the 
merger). In that calculation, the formation channels that do lead 
to CO+CO WD mergers are altogether different from our standard 
case. For the qceA = 0.125 model, binaries which end up merg- 
ing as CO+CO WDs generally start out with larger semi-latera recta 
(a x ( 1 .0 — e 2 )) of ~50 — 250 Rq , while for the standard model ini- 
tial semi-latera recta typically range from 20 — 100 Rq. For these 
merging WD pairs, the primary WDs are in general less massive 
than their standard model counterparts. 

The 0.2% of CO+CO WD merger progenitors which do en- 
counter the stage VIII accretion phase in the qceA = 0.125 model 
are born with near-equal masses, are closer at birth (semi-latera 
recta 20 — 100 Rq), and the first-formed WD always forms either 
during or after the CE (in contrast to our standard example, where it 
forms prior to the CE event). In all of these systems, there is a mass 
ratio reversal during the evolution, which is significant enough to 
allow for rapid rejuvenation of the secondary and the secondary is 
the first star to evolve into a WD (this process is also encountered 
on occasion in the standard model; Section [2~3l . However, these 
binaries do not meet our criteria for SN la progenitors: the most 
massive primary WD (at merger) in the 0.2% group is 0.7 Mq. 
This stems from the fact that the ZAMS progenitor masses for this 
channel are Quite l ow (~ 1.9 Mq). In any case, as demonstrated in 
IRuiter et all fcoilh the expected rate of merging CO+CO WDs is 
extremely low for this CE model. 

If we consider the '7— formalism' for CE evolution (e.g. 
iNelemans et~ai]l200d). which was assumed for all CE events for 
model G1.5 in lRuiter et al] ( l201ll) . stage VIII is also not encoun- 
tered very often. The typical formation channels which lead to dou- 
ble WDs from the standard model do not make double WDs in the 
'7— formalism' model, since for the latter the stars have larger or- 
bits post-CE, and will not merge within a Hubble time. Four % of 
primary CO WDs undergo stage VIII mass accretion, while for 3% 
it is the secondary WD for which the phase occurs. Both groups 
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Figure 8. As Fig. \5\ but using the WD masses as they were at birth. No 
matter which <j-cut is used, no theoretical Afbol distributions exhibit a 
shape similar to the absolute magnitude distribution of SNe la from lLi et alj 



have go = Jjl.2, and often the first interacting event is a CE, fol- 
lowed by a second CE after the primary has evolved into a CO WD. 
Both CE events enable the stars to be on a close enough orbit such 
that RLOF between a slightly evolved helium star and the primary 
WD is possible . In terms of f requen cy of events, the corresponding 
DDS model of iRuiter etal] < l201ll) yields a rate which is 10% of 
that of the standard model, thus we would expect a smaller fraction 
of potential SN la progenitors from violent mergers to be produced 
via the G1.5 model. 

Not surprisingly, binary population synthesis predictions for 
double degenerate mergers are rather sensitive to the assumed CE 
parametrization. In summary, a small reduction of CE efficiency 
(o;ceA = 1 — > 0.5) does not qualitatively affect our results. 
However a very low efficiency (oceA = 0.125) or use of the 7- 
prescription alone dramatically reduces both the total merger rate 
and the fraction of systems that undergo the stage VIII accretion 
phase. We note add itionally that using a higher CE-ejection ef- 
ficiency («ce > 1, iMeng & Yandr2012l) should lead to a larger 
number of WD mergers. 



6.2 WD mass growth: implications 

As an experiment to confirm that the mass-growth of the primary 
WD after it is born is key to our findings, in Fig. [8] we show re- 
sults equivalent to Figure[5]but using the StarTrack WD primary 
masses at their birth (dashed line in Figure [T) rather than at explo- 
sion. It is immediately obvious that the shape o f the observed ab- 
solute magnitude distribution of iLi et all d201 ll) is not reproduced 
in this case. Rather, for all q-cuts used, the theoretical Afboi distri- 
bution is rather flat from —17 to —19 mag, with a steep drop-off at 
A/boi = -19 mag. 

While Fig. [8] is only an experiment and not a self-consistent 
calculation for our StarTrack simulations, it may be roughly 
representative of predictions from alternative models. In particu- 
lar, we note that other population synthesis codes predict a total 
CO+CO WD merger mass distribution that is fairly smooth (Nele- 



mans 201(0; see also lBogomazov & Tutukovll20lbl fig. 1), in con- 
trast to the total mass distribution from StarTrack which has a 
second peak as a consequence of stage VIII (cf. Fryer et al. 2010 fig. 
3). We have identified an evolutionary channel that produces fairly 
massive primary WDs via helium-accretion, though it remains to be 
established whether such a prediction is confirmed by future obser- 
vations (see Footnote 3). In addition, improved input physics from 
models of interacting binaries that incorporate detailed stellar struc- 
ture - in particular for helium star models - are needed to properly 
add ress the issue of mas s transfer/accretion in close binaries (see 
e.g. lEggleton et alj|201 ll) . 



6.3 Ultra-prompt mergers 



As mentioned in section [5",2,l| the SNe la from WD mergers with 
delay times < 100 Myr ('ultra-prompt') are formed by way of a 
double CE: the initially less massive star loses its H-rich envelope 
while in the Hertzsprung gap or on the giant branch, and then later 
loses its He-rich envelope when it is a somewhat-evolved helium 
star. We find that 0.5% of all CO+CO mergers which have primary 
WD masses > O.9M0 merge less than 1000 yr after the last CE. 
The fraction is raised to 4% if the merger occurs within 10,000 
yr after the last CE. This has very interesting implications for at 
least a couple of reasons. The post-explosion spectra of these merg- 
ers are likely to show signatures of circumstellar material - thus 
far tho ught to be a likely sig nature of a single degenerate scenario 
SN la dSternberg et alEoi l|). On the other hand, since the primary 
WD will still be quite hot upon merging, these binaries may meet 
the necessary criteri a for producing 'core degenerate' mergers (see 
lllkov & Sokerll2012l) . which might produce SNe la - or some other 
bright transient - at a delay time larger than the merger time shown 
in Fig. [7] In any case, these types of mergers might contribute a 
small fraction to the population of sub-luminous SNe la. 



6.4 Violent vs. non-violent mergers 

IChen et all d2012l) used population synthesis to calculate rates of 
SNe la from CO+CO WD mergers and imposed some additional 
constraints, namely that in a dynamical merger some fraction of 
the total mass would be lost from the binary, in turn lowering the 
probability for a thermonuclear explosion. They found a possible 
upper limit of SNe coming from "violent mergers" to be about a 
factor of 5 lower than previous estimates from population synthe- 
sis. We woul d like to point out that by our definition, the mergers 
considered in lChen et alj d20K ) are n ot violent merger s. In the vi- 
olent mergers of IPakmor etaiT ( l2012l) - on which the IChen et alj 
d2012l) dynamically-motivated restrictions are based - no mass is 
lost in the merger itself. In the violent merger, the mass trans- 
fer phase is very short-lived, and the explosion takes place at the 
time of the merger at which point all of the matter is swept up 
jPakmor et alj2012T) . By contrast, the merge rs assumed to undergo 
mass-loss considered by IChen et all d2012h are likely to exhibit 
properties of non-violent mergers, where a hot envelope is accreted 
onto the primary WD. In such a case, a non-ne gligible amount of 
material mi ght be ejected f rom the binary (e.g. iFrver et al]|2010l 
but see also lDan et al,ll2012T) . though as already discussed, it is not 
clear that such systems will produce SNe la. 



8 Workshop on "Observational Signatures of Type la Supernova Progeni- 
tors", Leiden 2010 (http://www.lorentzcenter.nl/lc/web/ 1. 
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7 SUMMARY 

We find that the distribution of merging CO WD pairs predicted by 
the StarTrack binary evolution population synthesis code gives 
rise to a range of explosion brig htnesses via the prompt d etona- 
tion (violent merger) mechanism dPakmor et al .120 id. 1 20 lib that is 
compatible with observed SNe la (Fig.O- F urther, the DTD (Fig. [7) 
agrees well with the observed DTD fit from lGraur & Maozl d2012h . 
and as shown by Fig. [6] (lower panel), even gives rise to brighter 
SNe la preferentially at short del ay times - a trend th at is corrobo- 
rated by SN la observations (e.g. lWang & Han 2 o"l2h . 

We have iden tified tha t the g ood agreement with the bright- 
ness distribution o f iLi et alj J201 l|) depends critically on a particu- 
lar evolutionary phase during which the first-formed WD accretes 
mass from a companion which is a slightly evolved naked helium 
st ar (stage VIII , Fig. 2 ). In the standard population synthesis model 
of iRuiter et all ( 1201 ll) . nearly half of the CO+CO WD merger pro- 
genitors pass through such an accretion phase. 

Alternatively, if the critical He-rich accretion phase is not 
readily realised in nature, then the predicted explosion brightness 
distribution is not a good match to observations (Fig.[8]l, likely in- 
dicating that a different explosion scenario must dominate and drive 
the underlying shape of the SN la observed brightness distribution. 
Consequently, our results suggest that detailed studies of helium ac- 
cretion in binary systems will be a vital step toward establishing the 
plausibility of the violent merger model for a significant fraction of 
SNe la. 
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APPENDIX A: 

As mentioned in Section[2] we assume that when a non-degenerate 
star accretes from its companion, f a = 0.5, i.e. 50% of the mass 
lost from the companion is acc reted while the other 5 0% is lost 
from the binary (see sect. 3.4 of lBelczvnski et alj|2008h . In nature, 
it is not clear how much mass may be lost/gained in such a config- 
uration, and the assumed value (and variability) of f a has an im- 
portant impact on the e volutionary outcome (see the case study of 
iMennekens et alj|2010h . More important than the fraction of mass 
that is accreted might be in what manner angular momentum is lost 
from the binary system during non-conservative mass transfer (see 
also Claeys et al., in preparation, for a comparison study of four 
population synthesis codes). 

For the curren t study, we used th e published results from the 
standard model of iRuiter et alj d201 lh for which f a = 0.5 was 
assumed. However, to obtain a general idea of how much our re- 
sults might change if a different / a value had been used, we per- 
formed additional smaller-scale StarTrack runs for / a = and 
/ a = 1. We find that the relative frequency of CO+CO WD merg- 
ers is highest for the fully non-conservative case: 1.2 : 1 : 0.4 for 
a / a fraction of 0, 0.5 and 1, respectively. For more massive WD 
mergers however, the trend changes. For CO+CO WD mergers with 
-Mtot > 1.4M0, the relative frequency is 0.8 : 1 : 0.3 and for merg- 
ers which meet the criteria r\ = 1.5, m p > 0.8 Mq, the relative 
frequency is 0.5 : 1 : 0.3. 

Thus, assuming fully conservative mass transfer onto non- 
degenerate accretors (/ a = 1) seems at face value to result in 
an overall decrease in the number of CO+CO WD mergers. Our 
test here shows that the overall rates for massive CO+CO merg- 
ers might decrease at most by a factor of 3, though for fully non- 
conservative mass transfer the overall CO+CO WD merger rate 
could increase by 20%. The trend found he re is somewhat different 
than that found bv lMennekens et alj J20101) . however: one must be 
cautious and accept that full evolutionary calculations (including 
other stages of binary evolution) are needed to properly compute 
any relative change in our results, and such a parameter study is 
beyond the scope of this paper. 



APPENDIX B: 

Here we give the necessary information to reconstruct the fit in 
Figure [4] The 2-component fit was constructed with the following 
equality: 

7i max 

M bo l = ^ fl nmwD (Bl) 
n = 

where coefficients a n for the two different mass regimes can be 
found in Table B 1 . 

APPENDIX C: 

Since on average the most luminous SNe la are often found among 
relatively young stellar populations, while less luminous SNe la are 
found among older stellar populations (Sect. |5T2~Tt , we can explore 
how the shape of the violent WD merger peak brightness distribu- 
tion changes as a function of delay time. 

We split our WD mergers (from Fig. [5j into two age groups: 
delay times less than 1 Gyr and delay times greater than 1 Gyr. The 
corresponding peak brightness distributions for these 'younger' and 
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Table Bl. Coefficients for the polynomial and 
linear fits in Figure|4] 

0.8 < m W D < 1.15 



n 







2403.1822 


1 


-10943.254 


2 


19821.934 


3 


-17966.776 


4 


8140.1093 


5 


-1473.6569 


1.15 < m WD < 1.4 


n 


On 





-16.827815 


1 


-2.0807692 



'older' SNe la are shown in Fig. lCll and Fig. lC2l respectively. From 
the Figures, it is clear that the shapes of our theoretical brightness 
distributions (coloured lines) do not change very much as a function 
of age. Additionally, using a different time cut does not change the 
results much either, since both bright and dim violent mergers occur 
at early and late delay times (see Sect. l5.2l l. 

We additionally compa re our theoretic al distributions to the 
two luminosity functions of lli et all d20 1 lh that are separated by 
Hubble type. Li et al. (2011) split their luminosity functions into 
two groups based on host galaxy morphology: SNe la that are found 
in elliptical galaxies and spiral galaxies of type Sa, and SNe la that 
are found in spiral galaxies of type Sb and irregular galaxies (see 
their fig. 5). While delay times for individual SNe are not avail- 
able for the Li et al. (2011) sample, it is reasonable to assume 
that the average delay times for SNe la found among the group 
with elliptical galaxies will be larger than the average delay times 
for SNe la found in the group which contains irregular galaxies. 
This is reflected in the slightly higher proportion of low-luminosity 
events to high-lu minosity event s in Fig. IC2I compared to Fig. IC1I 
(greyscales). The lLi et"al] d201 lb sample (both age groups) contains 
SNe with different spectral properties, and we have included all 
74 objects regardless of their spectral classification (normal, high- 
velocity, 91T-like, 91bg-like, and 02cx-like). 

We note finally that the comparison of the theoretical bright- 
ness distributions to the observational data is for illustrative pur- 
poses only, since an accurate comparison requires detail ed infor- 
mation on the observed delay time distribution f rom the Li et al.l 



J201 lh sample. Furthermore, splitting the entire iLi et al 



feOllI) 



sample as we have done increases statistical noise: some bins only 
contain one SN event (see fig. 5 of Li et al. (2011)). Despite 
these potential caveats, the overall distribution shapes (observed 
and theoretical) show similar trends, and compare relatively well. 
In any case, it is possible (even likely) that other formation chan- 
nels and/or explosion mechanisms are also contributing at some 
level to the SN la p opulation at early and/or late delay times (e.g. 



level to tne aiN la p ( 
iKromer et ai]|2012l) . 
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Figure CI. Greyscale shows the peak brigh tness dis t ributio n of SNe la from 
Sb spiral and irregular galaxies as shown in lLi et aflfoOllL fig. 5). Coloured 
lines show our theoretical peak brightness distributions for WD mergers 
with delay times less than 1 Gyr. Such a split enables all 'prompt' SNe la 
from violent mergers to be included in the distribution, while still including 
some SNe with intermediate delay times, which would be expected in an 
observ ed sample invo lving spiral hosts (the 'Sb spiral & irregular' distribu- 
tion of lLi et alj \20\\\) ). 
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Figure C2. Greyscale shows the peak brightness distribution of SNe la from 

Li et all 1201 ll fig. 5). Coloured 
lines show our theoretical peak brightness distributions for WD mergers 
with delay times greater than 1 Gyr. 



